Loop structure of the lowest Bloch band for a Bose-Einstein condensate 
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We investigate analytically and numerically Bloch waves for a Bose-Einstein condensate in a 
sinusoidal external potential. At low densities the dependence of the energy on the quasimomentum 
is similar to that for a single particle, but at densities greater than a critical one the lowest band 
becomes triple-valued near the boundary of the first Brillouin zone and develops the structure 
characteristic of the swallow-tail catastrophe. We comment on the experimental consequences of 
this behavior. 

PACS numbers; 03.75.Fi, 05.30.Jp, 67.40.Db 



The possibility of studying experimentally the proper- 
ties of Bose-Einstein condensates in periodic potentials 
created by optical lattices (see e. g. (|, |§]) has stimu- 
lated a number of theoretical investigations. These in- 
clude nonlinear tunneling phenomena H, Bloch os- 
cillations 1^, 1^, stability of Bloch waves 0, and oscil- 
lations of condensates and localized excitations In 
a two-level model of nonlinear Landau-Zener tunneling, 
evidence was found for a loop structure in the adiabatic 
energy levels [|| . Recently Bronski et al. discovered a 
remarkable exact solution of the Gross-Pitaevskii equa- 
tion for a condensate in a one-dimensional optical lattice 
for a quasimomentum corresponding to the boundary of 
the first Brillouin zone. The reason for an exact solu- 
tion being possible is that the spatially dependent parts 
of the lattice potential and the internal field due to in- 
terparticle interactions cancel. Subsequently, in Rcf. ||lo[| 
it was shown by numerical studies that the current of 
the condensate for quasimomenta in the vicinity of the 
zone boundary had a singular behavior, and it was spec- 
ulated that the chemical potential of the condensate be- 
came multivalued near the zone boundary. In this Letter 
we use analytical and numerical methods to show that 
the lowest band has the structure associated with the 
swallow-tail catastrophe , and we comment on impli- 
cations of this result for experiment. 

The potential V{x) produced by a one-dimensional 
standing light wave has the form 



F(x) = K, sin^^. 



(1) 



where d is the period of the lattice. We shall assume that 
the condensate may be described by the Gross-Pitaevskii 
functional for the energy density, which when averaged 
over a lattice period is given by 
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and the mean number density n is given by 
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Here ip is the wave function of the condensate, m is the 
particle mass, and the effective two-body interaction is 
given by Uq = inh^a/m, where a is the scattering length 
for two-body collisions. The condensate wave function 
satisfies the nonlinear Schrodinger equation 



2m dx"^ 



+ [Vix) + Uo\M^-iL]i' = 0. 



(4) 



The quantity fi is the chemical potential, and this de- 
pends both on the mean density of particles, and on the 
quasimomentum. We shall look for solutions where the 
mean density is independent of position, and therefore 
these must satisfy the condition for quasi-periodicity. 



tpix + d) = exp(ifco?) ipix) 



(5) 



for all X. The quantity hk is the quasimomentum of 
the condensate. Our objective is to find the energy per 
particle, the chemical potential, and the current density 
of nonlinear Bloch waves as a function of k and of n. The 
current density of the condensate is given by 



Zmi \ dx dx 



It is independent of x provided satisfies Eq. 
direct calculation one can demonstrate that 
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The quantity j /n is the mean particle velocity. 

Without loss of generality, we may choose tp{0) to be 
real. Since the periodic potential is even in x, one infers 
from the quasi-periodicity condition (||) that the real part 
of the condensate wave function is an even function of x 
and the imaginary part an odd function. 

The solution at the zone center (fc = 0) corresponds to 
a purely real and periodic wave function, 'ip{x+d) — ip{x). 
Also at the zone boundary (fc — ±7r/c?) there is a purely 
real solution which is antiperiodic, ip{x + d) = —ip{x). 
Since ip is real, these solutions have zero current, j = 
0. Other solutions are complex and the current is non- 
zero. Since the potential in Eq. (m is real, it follows that 
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'4'kix) = tp-ki^)- Naturally, all physical quantities are 
periodic in k with period 2-k /d. 

The structure of the energy as a function of k near the 
zone boundary may be brought out by making a varia- 
tional calculation of the energy. We use a trial function 
of the form 

'ip{x) — ^/n{cosa expi/cx + sino; exp(z[fc — 27r/d]a::)}, (8) 

where a is a variational parameter. This is similar to that 
used in the almost-free-electron approximation for the 
band structure of a single particle, and it also corresponds 
to the two- level model studied in Ref. Q . When the trial 
function (^) is inserted into the energy functional (||), it 
becomes 
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where the three terms are the kinetic, potential and inter- 
action energies, respectively, with I = {ck + efc-27r/d)/2, 



where — h'^k^/2m, and Ae 



{2TTh /md){k ~ n/d). The optimal value of a is obtained 
by requiring (||) to be stationary with respect to varia- 
tions in a, and the chemical potential /x is then found by 
differentiating the resulting energy with respect to the 
density, ^ — d£/dn. 

At the zone boundary k = ztn/d the stationarity of (||) 
yields two solutions: sin 2a = Vo/2nUo and cos 2a — 0. 
The first solution, which requires that the density exceed 
a critical value ric given by 



= Vo/2Uo, 



has an energy density 
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a chemical potential 
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where 

and a current density 
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The two possible signs for the current density reflect the 
fact that if a is a solution to the stationarity condition, 
so is 7r/2 — a. 



Equations (10)-(|14|) are in fact exacts the variational 
trial function TS) being in this case identical with the 



analytical solution found in Refs. H, |T^. At the critical 
density when the analytical solution first appears with 
a = 7r/4, the current is zero since il){x) is real. We note 
that the critical chemical potential is the height of the 
barrier, plus the kinetic energy of a free particle with 
a momentum equal to the quasimomentum at the zone 
boundary. 

For n > Tic the exact solution (H) with sin 2a = 
Vo/2nUo < 1 is complex and carries a non-zero current. 
Remarkably, although this state has k = ztir/d, it has 
neither the highest energy nor the highest chemical po- 
tential of all states in the band: the real antiperiodic 
wave function, also with k = ±11 /d, has higher yit, and 
it corresponds to the second solution of the stationarity 
condition, cos 2a = or a = 7r/4 (the other possibil- 
ity, a — 37r/4, belongs to the next, excited Bloch band 
and will not be considered here). This solution exists 
at all densities and corresponds to a chemical potential 
fJ'2 = fJ-c + 3(n — nc)Uo/2 which is higher than fii at 
n > Tic- 

For a general value of fc, the Ae term in Eq. (|^) 
for the total energy tends to make the roots of the 
stationarity condition initially at a — 7r/4 and a = 
[arcsin(Vo/2n?7o)]/2 approach one another, and the two 
roots merge and disappear for a sufRciently large value of 
the magnitude of Ae. The energy per particle obtained 
using the trial function (^) is shown in Fig. 1 , where one 
sees the swallow-tail structure near the zone boundary. 
Throughout, we shall measure energies in units of huj, 
where uj = {2Va/m)^/'^i: / d is the angular frequency of 
small oscillations about a minimum of the potential. 
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FIG. 1: Swallow-tail structure of the energy per particle as a 
function of k for Vb = ^fiuj and n/nc= 1.20. 

A quantitative analysis shows that the trial form (||) is 
a good approximation near the zone boundary provided 
h'^'K^ /2md^ ^ nf/o- The analytical calculations may be 
improved by including plane- wave components with wave 
numbers that differ from the ones in Eq. (|^) by multiples 
of the reciprocal lattice vector 2-K/d. 

We now demonstrate that the swallow-tail structure 
also arises in numerical solutions of the full Gross- 
Pitaevskii equation (El). At condensate densities below 
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the critical one, the wave function and dispersion rela- 
tion for nonlinear Bloch waves are qualitatively similar 
to those for a single particle. To illustrate this, we again 
consider the case where the potential barrier is Vo — 3huj. 
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FIG. 4: The current density j as a function of quasimo- 
mentum at densities n/nc—0.57 (dotted), 0.78 (dashed), and 
0.99 (full line). Here jo — nhn/md is the current density at 
k — Tr/d in the absence of a lattice. 



FIG. 2: Chemical potential /i in units of hu; as a function 
of quasimomentum for densities n/ric— 0.57 (dotted), 0.78 
(dashed), and 0.99 (full line). 



analytical solution. 



Vo nx 
Uq d 



(15) 
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FIG. 3: Energy per particle, £/n, in units of hui as a function 
of quasimomentum for densities n/nc= 0.57 (dotted), 0.78 
(dashed), and 0.99 (full line). 



The current density j(fc, Uc) goes to zero at A: = zLir/d but 
has an infinite derivative at this point, see Fig. 4. Conse- 
quently, ^,{k,nc) develops a singularity at A: = ivr/d (see 
Fig. 2), as does the energy density. 
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In Figs. 2 and 3 we plot the chemical potential /i(fc, n) 
and the energy per particle for several values of density 
n < Uc- In all cases the chemical potential lies below the 
top of the potential barriers. Therefore, the energy bands 
are relatively narrow as the motion of the condensate 
involves tunneling through the barriers. The difference 
between the chemical potential and the energy per par- 
ticle reflects the nonlinear character of the interparticle 
interaction. 

In Fig. 4 we show the current density j{k, n) of the non- 
linear waves for the same number densities as in Figs. 2 
and 3. The current goes to zero both at the zone cen- 
ter and the zone boundary. At the critical density the 
real antiperiodic wave function corresponding to the top 
of the band and having k — ±n/d coincides with the 



FIG. 5: Chemical potential as a function of quasimomentum 
for densities n/ric — 1.14 (dashed) and 1.25 (full line). 

For n > Uc the Bloch wave evolves from a real peri- 
odic wave function at fc = to the analytical plane-wave 
solution (H) at k = n/d, which, however, is now com- 
plex. At k = n/d there is a second solution equal to 
the complex conjugate of the first one and, in addition, 
a third one which is real, and which therefore carries no 
current. Also for k close to Tr/d one finds three solutions 
for a given k. We plot the chemical potential fi{k,n), 
the current density j(k,n)^ and the energy per particle 
£(fc, n)/n in Figs. 5-7, respectively, for two values of the 
density n > ric- We note that j(fc, n) calculated from the 
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wave function using Eq. (m is consistent with the relation 
(I). 
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FIG. 6: Current density j as a function of quasimomentum 
for densities, n/ric = 1-14 (dashed) and 1.25 (full line). The 
unit of current jo is the same as in Fig. 4. 
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FIG. 7: Energy per particle f /n of a Bloch wave as a function 
of quasimomentum for a condensate density n/uc = 1.20. 

The unusual structure of the energy bands has impor- 
tant implications for experiment. If a single particle in a 
periodic potential is subjected to a small constant force, 



k varies linearly in time, and the particle velocity and po- 
sition execute Bloch oscillations. Such behavior has been 
observed in optical lattices for thermal atoms above the 
Bose-Einstein transition temperature Similar argu- 
ments apply to the motion of a Bose-Einstein conden- 
sate, provided the energy of the condensate in the lowest 
band is a continuous, single-valued function of k. Bloch 
oscillations have indeed been observed for condensates 
Q. At densities greater than ric, the situation is differ- 
ent because of the swallow-tail structure of the energy 
band. When a weak force is applied to a condensate ini- 
tially at rest, in the simple picture of the motion k will 
increase and with time will reach the value Tr/d. With 
further increase of fc, the system will continue on the 
branch of the spectrum for k > n/d which has the same 
slope at fc = Tr/d as that of the initial branch. With this 
choice of branch, the wave function varies continuously 
as k increases past Tr/d, whereas there is a discontinuous 
change in the wave function of the lowest energy state. 
A dramatic manifestation of this is that the current den- 
sity in the lowest energy state changes sign at A; = n/d. 
With time, k will eventually reach the tip of the swal- 
low tail, and it then becomes impossible to describe the 
condensate motion in terms of the usual adiabatic pic- 
ture, as was found in Ref. ||^ for a two-level model. How 
to treat the dynamics under such conditions remains a 
challenging open problem. 

The experiments on Bloch oscillations in Ref. ||] were 
carried out in a regime in which the interaction energy 
was small compared with the depth of the potential due 
to the optical lattice, and therefore one would not have 
expected the energy band to have the swallow-tail struc- 
ture. It is important to extend such experiments into the 
regime in which the strong nonlinear effects predicted 
here should occur. Theoretically, the stability of a con- 
densate moving in a periodic potential has been examined 
for relatively low values of the two-body interaction , 
and such studies should be extended to higher values of 
the interaction, where the swallow tail develops. 

We are grateful to D. Aristov and S.J. Chang for very 
useful discussions. 
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